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13  ABSTRACT  (Maximum  200  wonts) 

Presented  here  is  the  result  of  a  Finite  Element  Analysis  (FEA)  and  its  correlation  with  test  data  for  a  Class  IV  fiexten- 
sional  underwater  acoustic  transducer.  The  thick,  elliptical  fiberglass/epoxy  shell  of  the  transducer  Finite  Element  Model 
(FEM)  was  modeled  with  20-node  solid  elements  using  effective  engineering  properties  derived  from  a  micromechanic.,  theory. 
Unidirectional  gap  elements  were  used  to  model  the  clearance  between  the  end  curves  of  the  shell  and  the  two  D-shapcd  alumi¬ 
num  inserts.  The  D-inserts  were  placed  along  the  mayor  axis  of  the  shell  to  transmit  the  driving  power  of  the  piezoelectic 
ceramic  stack  in  order  to  excite  shell’s  flaxural  modes  for  an  efficient  acoustic  radiation.  The  model  also  featured  large  defor¬ 
mation.  The  predicted  global  displacements  of  the  shell,  as  well  as  the  local  strain  values  at  several  critical  locations  in  the 
shell,  correlated  very  well  with  the  test  data  obtained  from  the  stack  preload  and  hydrostatic  pressure  tests.  Gapping  between 
the  end  curve  of  the  shell  and  the  D-insert,  which  was  detected  during  the  pressure  test,  was  also  predicted  in  the  analysis. 
Gapping  has  been  shown  to  be  detrimental  to  the  structural,  as  well  as  the  acoustic  performance  of  the  shell. 

The  success  of  this  pilot  model,  in  accurately  predicting  the  structural  behavior  of  the  transducer  shell  under  simulated 
loading  conditions,  has  lead  to  a  series  of  FEA  studies  to  determine  the  sensitivities  of  the  shell’s  structural  performance  to 
thickness  variation  and  to  nonuniform  material  property  distribution  in  the  shell.  Nonuniform  material  property  distribution 
arise  from  the  variation  in  parameters,  such  as  moduli  of  the  constituents,  fiber  volume  fraction,  and  ply  arrangement  of  the 
laminates.  The  model  is  currently  being  used  to  investigate  the  concept  of  flexible  D-insert  designs  conforming  to  the  deformed 
shape  of  the  shell  in  order  to  maintain  a  closed  shell-to-D-insert  contact  for  deeper  underwater  applications. 
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ABSTRACT 

Presented  here  Is  the  result  of  a  Finite  Element  Analysis  (FEA)  and  Its  correlation  with  test  data 
for  a  Class  IV  flextenslonal  underwater  acoustic  transducer.  The  thick,  elliptical  fiberglass /epoxy 
shell  of  the  transducer  Finite  Element  Model  (FEM)  was  modeled  with  20-node  solid  elements 
using  effective  engineering  properties  derived  from  a  micromechanlcs  theory.  Unidirectional  gap 
elements  were  used  to  model  the  clearance  between  the  end  curves  of  the  shell  and  the  two 
D-shapcd  aluminum  inserts.  The  D-inscrts  were  placed  along  the  major  axis  of  the  shell  to 
transmit  the  driving  power  of  the  piezoelectric  ceramic  stack  in  order  to  excite  shell's  flexural 
modes  for  an  efficient  acoustic  radiation.  The  model  also  featured  large  deformation.  The  predicted 
global  displacements  of  the  shell,  as  well  as  the  local  strain  values  at  several  critical  locations  in 
the  shell,  correlated  very  well  with  the  test  data  obtained  from  the  stack  preload  and  hydrostatic 
pressure  tests.  Gapping  between  the  end  curve  of  the  shell  and  the  D-lnsert,  which  was  detected 
during  the  pressure  test,  was  also  predicted  in  the  analysis.  Gapping  has  been  shown  to  be 
detrimental  to  the  structural,  as  well  as  the  acoustic  performance  of  the  shell. 

The  success  of  this  pilot  model,  in  accurately  predicting  the  structural  behavior  of  the 
transducer  shell  under  simulated  loading  conditions,  has  lead  to  a  series  of  FEA  studies  to 
determine  the  sensitivities  of  the  shell's  structural  performance  to  thickness  variation  and  to 
nonuniform  material  property  distribution  in  the  shell.  Nonuniform  material  property  distribution 
arises  from  the  variation  in  parameters,  such  as  moduli  of  the  constituents,  fiber  volume  fraction, 
and  ply  arrangement  of  the  laminates.  The  model  is  currently  being  used  to  investigate  the 
concept  of  flexible  D-lnsert  designs  conforming  to  the  deformed  shape  of  the  shell  in  order  to 
maintain  a  closed  shell-to-D-lnsert  contact  for  deeper  underwater  applications. 

INTRODUCTION 

For  the  past  few  years,  the  Research,  Development,  Test  and  Evaluation  Division  of  Naval 
Command,  Control  and  Ocean  Surveillance  Center  (NRaD),  has  been  tasked  to  develop  a  high- 
power,  low-frequency  transducer  as  an  effective  underwater  sound  source.  A  2-lnch  thick 
fiberglass /epoxy  shell  with  elliptical  cross  section  has  been  developed  recently  for  the  projector 
element  of  the  Class  IV  flextenslonal  design.  In  developing  this  transducer  to  achieve  adequate 
strength  in  the  shell  and  durability  in  service,  it  is  imperative  to  test  the  strength  of  the  composite 
shell  under  simulated  loading  conditions.  In  light  of  the  complicated  process  of  manufacturing 
composite  shells,  it  is  virtually  impossible  to  bring  important  parameters  such  as  shell  thickness 
and  the  material  strength  distribution  within  the  shell  to  typically  allowable  tolerances  for  the 
metallic  parts.  Due  to  the  random  nature  of  these  variations  and  high  cost  of  manufacturing  the 
composite  shells,  it  is  impractical  to  produce  enough  specimens  for  strength  tests  across  the  full 
spectrum  of  these  variables.  Therefore,  an  analytical  approach  is  necessary  to  complement  a  few 
control  tests. 

In  many  FEA  programs  that  can  analyze  the  composite  structures,  the  Classical  Lamination 
Theory  (CLT)  is  used  to  "smear"  the  ply  material  properties  of  composites  so  that  the  FEA  program 
secs  a  homogeneous  body  with  equivalent  orthotropic  properties.  Programs  that  use  CLT  can 
calculate  ply-by-ply  results,  but  they  cannot  take  into  account  through -thickness  effects  (in-planc 


normal  and  Interlaminar  shear  stresses).  MSC/NASTRAN  ( 1989)  developed  a  single  plate  or  shell 
element  with  through-thickness  shear  chert  by  incorporating  the  equilibrium  conditions  of  a  'bent 
txam".  However,  this  shell  element  is  suitable  for  thin  shell  application,  hut  not  for  the  thick  shell 
such  as  the  Class  IV  flextenslonal  shell  A  study  on  limitations  and  shortcomings  of  using 
MSC/NASTRAN’s  composite  shell  elements  for  analysis  of  thick  shells  was  reported  by  Shaw 
11989).  Using  three-dimensional  (3-D)  solid  elements  is  an  alternative  to  model  the  composite 
structure.  Rouse  (1988)  pointed  out  that  3-D  solid  elements  coupled  with  lamination  theory  to 
model  laminated  structures  has  become  more  popular  among  the  composite  designers.  The  use  of 
"smeared"  material  properties  for  the  solid  elements  means  fewer  elements  are  needed  to  model 
the  structure,  but  ply  data  are  still  available  at  the  end  of  analysis.  This  approach  is  chosen  here 
to  analyze  flextenslonal  transducer  shells.  The  ABAQUS  FEA  program  is  executed  on  a  Convex 
mini-supercomputer.  The  analysis  results  are  validated  through  comparisons  with  test  data 
obtained  from  stack  preload  and  hydrostatic  pressure  tests.  The  data  compared  include  the  major 
and  minor  axis  deflections  and  strain  values  at  several  critical  locations  in  the  shell.  The  natural 
frequencies  predicted  by  the  model  are  also  compared  with  the  frequency  data  obtained  from  the 
modal  survey  on  the  bare  shell. 

FINITE  ELEMENT  ANALYSIS  OF  COMPOSITE  SHELL 


MODEL  GEOMETRY 

Two  5.7-inch  long,  2-lnch  thick  (nominal)  fiberglass /epoxy  shells  were  fabricated  for  the  tests. 
The  test  configuration  for  the  test  specimen  (Shell  010)  showing  strain  gage  locations  is  giving  in 
Figure  1.  In  constructing  the  FEM  of  the  composite  shell,  the  inner  surface  of  the  shell,  which 
consists  of  two  radii  (at  "flat"  and  "curve")  meeting  at  a  transition  point,  was  chosen  as  the 
reference  surface  on  which  the  geometry  of  the  model  was  defined.  Two  cylindrical  coordinate 
systems  originated  at  the  centers  of  the  "flat"  and  "curve"  were  used  to  facilitate  the  generation  of 
the  shell  model  using  the  PATRAN  FEM  pre  and  postprocessors  hosted  on  a  MlcroVAX  1I/GPX 
workstation.  Figure  2  shows  a  1/8-symmetric  (1 -quadrant  and  1  /2-lcngth),  3-D  FEM  of  the  test 
specimen  constructed  with  ABAQUS’  20-nodc  solid  elements  with  2x2x2  reduced  integration 
points  (C3D20R).  The  shell  is  Interfaced  at  the  inner  curve  by  a  4 -inch  long,  D-shapcd  aluminum 
insert  which  was  also  modeled  with  C3D20R  solid  elements.  The  D-lnscrt  has  a  matching  radius 
slightly  smaller  than  the  end  radius  of  the  shell.  Thus,  a  clearance  exists  along  the  D-lnscrt  and 
shell  interface,  as  shown  in  Figure  3.  This  initial  clearance  was  modeled  with  ABAQUS 
unidirectional  gap  elements  (GAPUNI)  defined  by  a  set  of  matching  pairs  of  D-lnscrt  and  shell 
nodes  along  the  radial  direction  of  the  D-lnsert.  The  presence  of  these  gap  elements,  which  close 
and  transmit  compressive  forces  in  the  radial  direction  of  the  D-inscrt  when  loaded,  makes  the 
shell  model  a  FEM  with  the  "boundary  nonlinearity".  Note  that  current  analyses  also  consider  tnc 
"geometric  nonlinearity"  to  account  for  the  possible  large  deformations  in  the  loaded  shell. 

"SMEARED"  MATERIAL  PROPERTIES 

The  laminated  construction  of  the  shell  wall  in  the  current  Class  IV  transducer  design  can  be 
considered  as  a  balanced,  symmetric,  angle-ply  laminate  (l.c.,  there  arc  equal  numbers  of  plies 
with  an  equal  but  opposite  orientation  angle  of  7C).  Instead  of  performing  a  time-consuming, 
ply-by-ply  material  and  structural  definition  for  the  composite,  the  effective  engineering  properties 
of  the  laminate  were  estimated  from  the  referenced  micromechanical  equations.  These  equations, 
expressing  the  ply  properties  in  terms  ol  their  constituent  properties  and  volume  fractions,  arc  ( 1 ) 
the  rule  of  mixtures  (sec  Halpin,  1984)  to  estimate  the  longitudinal  modulus  (Ej  j)  and  Poisson's 
ratio  (vj2),  and  (2)  the  Halpin-Tsai  empirical  formulas  (see  Halpin.  1984)  to  approximate  the 
transverse  and  ln-planc  shear  moduli  (E22  and  G12  respectively).  To  account  for  the  variation  in 
the  fiber  volume  fraction  through  the  thickness,  the  values  of  "smeared"  orthotropic  material 
properties  for  the  test  specimen  were  evaluated  for  a  range  of  fiber  volume  fractions  between  55 
and  65  percent.  Best  estimates  of  the  upper-  and  lower-bound  values  based  on  results  of  that 
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evaluation  were  provided  by  Maltby  (1988)  and  are  tabulated  in  Table  1  The  material  orihotmpy 
was  accounted  for  by  orienting  the  material  properties  of  the  elements  in  the  "Oat"  and  'curve 
portions  of  the  shell  to  the  respective  local  cylindrical  coordinate  system  defined  for  the  model,  a-, 
shown  in  the  figure  below  Table  1.  However,  only  linear  and  homogeneous  proper!  ies  were  used  m 
the  analyses. 

BOUNDARY  AND  LOADING  CONDITIONS 

The  symmetry  conditions  In  X-Y,  Y-Z  and  Z-X  planes  shown  in  Figure  2  provide  adequate 
constraints  on  the  model  in  the  modal  analysis  to  calculate  natural  frequencies  of  the  shell  and  in 
the  stress  analysis  to  predict  the  response  of  the  shell  due  to  the  stack  preload  and  the  hydrostaic 
pressure.  Before  pressure  is  applied  in  the  hydrostatic  pressure  test,  the  shell  Is  preloaded  in  the 
stack  to  136  kips.  To  simulate  this  loading  condition,  the  analysis  is  performed  in  two  steps.  In 
step  1,  a  prescribed  displacement  In  the  outward  direction  of  the  major  axis  (X-direction)  is 
imposed  on  the  wedge  end  of  the  extension  rod  (node  2073  in  Figure  2)  such  that  the  reacting 
force  In  the  rod,  which  replaces  the  ceramic  stack  in  the  test  configuration  and  is  modeled  as  a 
spring  element  with  a  stiffness  equivalent  to  that  of  the  extension  bar  and  the  wedge  assemiy  in 
series  (K  =  2. 1438  x  10®  lb/In).  reaches  the  value  of  136  kips.  Note  that  a  constraint  equation  in 
X-dlrection  is  Imposed  among  the  nodes  at  the  slot  of  the  D-insert  (nodes  1639,  1647.  1655, 

1900,  1963,  2026,  etc.)  to  ensure  that  the  constraint  force  induced  at  node  2073  due  to  the 
precribed  displacement  Is  properly  transmitted  to  the  D-insert.  Also  note  that  several  Iterations 
are  necessary  to  reach  the  desirable  value  of  the  preload  in  step  1 .  In  step  2.  the  pressure  is 
applied  at  the  external  surface  of  the  shell  while  this  prescribed  displacment  is  maintained,  even 
though  the  curve  end  of  the  shell  deflects  outward  under  increasing  pressure.  This  results  in  a 
relaxation  of  the  stack  load. 

Table  1.  Estimated  Engineering  Properties  for  Shell  010 


’’Smeared' 

'  Engineering  Properties 
(xTo®psi) 

Poisson's 

Ratios 

En 

E22 

e33 

g12 

g13 

g23 

v12-v13- v23 

Lower- Bound 

5.7 

1.8 

1.2 

0.57 

0.57 

0.57 

0.3 

Upper-Bound 

7.0 

2.4 

2.3 

0.90 

0.90 

0.90 

0.3 

Coordinate  Systems  For 
Engineering  Properties 


Notes:  (1)  The  engineering  constants  defined  above  are  transformed  into  R  (radial). 

©(circumferential),  and  Z  (axial)  coordinate  system  used  in  ABAQUS  analyses. 
The  indices  1,  2,  and  3  in  Table  1  become  2.  3.  and  1  in  ABAQUS. 

(2)  The  conditions  relating  Poisson’s  ratios  to  Young's  moduli,  Vy/Ey  =  v^/Ey 
(i,  j  =  1 ,  2,  and  3),  apply  when  reversing  indices  of  vy 
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ANALYSIS  RESULTS  AND  COMPARISON  WITH  TEST  DATA 


MODAL  ANALYSIS 

The  calculated  natural  frequencies  of  the  two  models  encompass  the  test  result  obtained  from 
a  modal  survey  of  the  bare  shell,  with  the  lower-bound  property  model  rendering  better 
correlation.  The  comparisons  are  given  in  Table  2.  The  lower-bound  property  model  predicts  the 
lowest  natural  frequency  to  be  at  164.41  Hz,  which  Is  within  3.9  percent  of  the  test  result  of  17] 
Hz  obtained  in  a  modal  survey  by  Maltby  (1991).  The  upper-bound  property  model  overpredicts 
the  lowest  natural  frequency  by  8.4  percent  at  185.33  Hz.  The  predicted  mode  shapes  for  the  first 
two  modes  are  shown  In  Figures  4a  and  4b.  No  mode  shape  data  from  the  modal  survey  was 
available  for  comparison. 


Table  2.  Dynamic  Characteristics  of  Shell  010, 
Modal  Analysis  Versus  Modal  Survey  Data 


Natural  Frequencies  (Hz) 


1st  Mode 

2nd  Mode 

3rd  Mode 

Hz  Test  Comp 

Hz 

Hz 

Lower-Bound  Model 

164.41  -3.9% 

840.50 

1,757.6 

Upper-Bound  Model 

185.33  +8.4% 

960.29 

1.950.5 

Modal  Survey  Data 

171. 

As  shown  in  Figure  4a,  the  opposite-phased  bending  of  the  "flat"  and  the  "curve”  characterizes 
the  lowest  mode  of  free  vibration  for  the  oval  shell.  i.e„  the  "flat"  bends  inward  and  the  "curve" 
displaces  outward.  Although  the  concept  for  flextenslonal  transducers  is  based  on  the  mutual 
coupling  between  two  perpendicular,  flexural  modes  of  vibration  for  an  oval  shell,  the  Interaction 
between  the  membrane  and  the  flexural  modes  plays  an  important  role  in  determining  the 
structural  behaviors  of  the  shell. 

STACK  PRELOAD  TEST  RESULTS 

The  prescribed  displacements  at  node  2073  to  Induce  a  136-kip  stack  preload  in  the  first  load 
step  are  0.2654  inch  and  0.2330  inch  for  the  lower-  and  upper-bound  property  models 
respectively.  The  stack  loads  are  plotted  against  the  major  and  minor  axis  displacements  and  the 
results  are  compared  with  the  test  data,  as  shown  in  Figures  5a  and  5b.  As  can  be  seen,  excellent 
correlations  are  achieved,  as  the  test  results  fall  between  the  predicted  values  of  the  two  models. 

HYDROSTATIC  PRESSURE  TEST  RESULTS 

The  prescribed  displacement  derived  In  the  first  load  step  Is  maintained  when  external  pressure 
is  applied  in  the  second  load  step.  The  flattening  of  the  shell  under  external  pressure  relaxes  the 
stack  load  in  the  D-lnsert.  The  calculated  results  of  pressure  load  versus  stack  load  are  plotted 
against  test  data  and  shown  in  Figure  6.  Fairly  good  correlation  between  the  calculated  and  test 
results  is  seen.  Note  that  each  calculated  data  point  In  Figure  6  represents  a  convergent  solution 
of  the  Incremental  loading  procedure  carried  out  In  ABAQUS.  A  total  pressure  of  459. 1  psi  and 
four  Increments  are  specified  at  the  beginning  of  the  second  load  step  In  the  analysis.  As  planned, 
convergent  solutions  are  obtained  In  four  equal  Increments  for  the  upper  bound  model.  However, 
after  convergent  solutions  are  obtained  at  the  first  two  equal  Increments  for  the  lower-bound 
property  model,  the  size  of  the  pressure  increment  was  reduced  by  one-half  automatically  to  reach 
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a  solution  at  a  higher  pressure  level  in  ABAQUS  computer  run  In  all,  six  load  increments  an 
necessary  to  reach  a  solution  for  459. 1  -psi  pressure.  As  shown  in  Figure  6,  an  almost  identical 
linear  pressure- to-stack  load  relationship  exists  for  the  lowei-  and  upper -bom id  property  models, 
up  to  three-fourth  of  the  way  (344.3  psi)  of  the  maximum  pressure.  But  the  relationships  become 
nonlinear  and  diverge  for  the  two  models  when  pressurized  beyond  344.3  psi.  Tills  observation 
indicates  that  some  kind  of  nonlinear  behavior,  other  than  the  nonlinear  material  behavior,  must 
have  happened  during  the  last  increment  of  pressure  loading  that  changed  the  load  paths  and 
caused  load  redistribution  in  the  shell  models.  This  question  will  be  answered  later. 

The  minor  axis  displacement  versus  pressure  load  curve  is  plotted  in  Figure  7.  It  seems  that  the 
analysis  models  underestimate  the  flexural  stiffness  of  the  shell.  The  minor  axis  deflection  for  the 
upper-bound  property  model  is  0.289  inch  at  the  maximum  pressure  of  459. 1  psi.  as  compared  to 
0.176  inch  for  the  test,  a  64  percent  discrepancy.  Three  important  observations  are  maue 
regarding  results  shown  in  Figure  7: 

1.  The  unloading  curve  does  not  follow  the  loading  curve  for  the  test.  This  indicates  that  some 
nonlinear  material  response  has  taken  place  during  the  pressure  test.  It  should  be  noted  that 
the  unloading  curve  obtained  in  the  analytical  simulation  of  unloading  follows  the  loading 
curve  almost  exactly.  This  result  is  not  surprising  since  linear-elastic  material  properties 
were  used  in  both  analysis  models. 

2.  The  pressure  load  and  displacement  relationship  becomes  nonlinear  when  the  shell  is 
pressurized  beyond  344.3  psi.  This  nonlinear  behavior  is  the  result  of  permitting  large 
deformations  in  the  current  analysis. 

3.  The  flexural  behavior  is  sensitive  to  the  variations  of  the  material  properties,  especially  when 
the  shell  is  subjected  to  high  pressure. 

GAPPING  OF  SHELL 

A  closer  look  at  the  progression  of  deformation  resulting  from  incremental  loading  of  pressure 
icveals  that  gapping  has  occurred  at  the  shell  and  D-lnsert  interface  at  the  last  load  increment 
The  shell  and  D-insert  interface  remains  closed  during  the  first  three  increments  of  pressure 
loading  (up  to  344.3  psi),  but  opens  and  forms  a  gap  near  the  end  of  the  shell  when  the  external 
pressure  reaches  459. 1  psi,  as  shown  in  Figure  8.  An  enlarged  view  clearly  depicting  this  interface 
gapping  is  giving  in  Figure  9.  It  is  important  to  report  that  this  behavior  was  also  observed  in  the 
pressure  test.  Gapping  can  be  attributable  to  an  Interface  condition  that  develops  between  the 
shell  and  the  D-insert.  As  pressure  load  increases  and  flattens  the  shell,  the  end  curve  of  the  shell 
takes  on  a  parabolic  shape  and  its  radius  of  curvature  eventually  becomes  smaller  than  the  D 
insert  and  separation  occurs  at  the  apex  of  the  curve.  This  observation  signals  the  presence  of 
large  deformation  at  the  curve  of  shell.  Gapping  starts  from  the  apex  of  the  end  curve  and 
propagates  upward  along  the  interface  and  would,  eventually,  cause  the  shell  to  pivot  about  the 
upper  edge  of  the  D-lnsert  if  higher  pressure  is  applied.  This  phenomenon  is  undesirable  for  three 
reasons  as  follows: 

1.  The  concept  of  a  flextensional  transducer  relies  on  the  premise  that  a  close  coupling  exists 
between  the  ceramic  stack  length  mode  and  shell  flexural  modes.  In  the  most  general  sense, 
the  stack  extends,  transmitting  its  load  through  the  D-inscrts  to  the  shell,  which  in  turn 
flexes,  thus  driving  the  name  flextensional.  The  movement  of  the  contact  area  away  from  the 
apex  of  the  shell  moves  the  path  of  the  stack  load  away  from  the  major  axis  of  the  shell  and 
reduces  the  effectiveness  of  the  D-lnscrt  in  transmitting  the  maximum  load  from  the  stack  to 
the  shell.  The  reduction  of  load  tramsmitted  by  the  D-inscrt  undoubtly  decreases  the  bending 
moment  generated  to  drive  the  flexural  modes  of  the  shell  vigorously  and.  therefore,  can 
reduce  the  acoustic  performance  of  the  shell. 
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2.  Gapping  causes  redistribution  of  contact  load  along  the  interface  to  a  less  favorable 
position.  The  contact  stresses  are  more  concentrated  on  the  area  of  shell  that  maintains 
contact  with  the  upper  edge  of  D-insert.  The  significant  stress  components  for  the  last  load 
increment  of  the  upper-bound  property  model  are  given  in  Figures  10  to  13.  The  stress 
component  most  affected  by  the  change  of  Interface  conditon  is  the  radial  stress  (S33).  As 
shown  in  Figure  10.  the  gradients  of  this  stress  distribution  becomes  more  abrupt  when  the 
interface  condition  changes.  The  maximum  radial  stress  is  10,147  psi  in  compression,  which 
occurs  at  the  narrow  strip  of  shell's  inner  surface  area  that  remains  in  contact  with  D-insert. 
However,  the  areas  immediately  surrounding  the  contact  are  in  tension  (up  to  2,386  psi)  due 
to  the  discontinuity  of  contact.  Although  the  magnitude  of  this  through-thickness  normal 
tensile  stress  is  not  very  high,  it  is  detrimental  considering  the  long  fatigue  life  demanded  of 
the  shell. 

3.  The  oval  shell,  which  is  supported  firmly  at  its  curve  ends  by  the  D-inserts,  is  tuned  to 
vibrate  at  a  certain  frequency  band  during  the  operation  of  a  transducer.  But,  gapping 
between  the  shell  and  the  D-insert  caused  by  high  depth  pressures  alters  the  end-support 
condition  of  the  shell.  As  shown  in  Figure  9,  the  gap  widens  and  the  shell  is  supported  by 
approximately  the  upper  one-fifth  of  the  D-lnsert  when  the  pressure  reaches  459. 1  psi.  The 
reduction  of  the  end-supporting  area  for  the  shell  changes  the  dynamic  characteristics  of 
the  shell  and,  therefore,  affects  the  resonance  condition  of  the  shell  during  operation  of 
the  transducer  at  greater  depths. 

HOOP  AND  SHEAR  STRESS  DISTRIBUTIONS 

Other  significant  stress  components  whose  distributions  are  not  affected  much  by  the  presence 
of  the  gap  are  the  hoop  stress  (S  j  and  the  through-thickness  shear  stress  (S13).  Figure  1 1 
shows  the  distributions  of  the  hoop  stress  for  the  maximum  pressure  of  459. 1  psi.  The  hoop 
stresses,  which  are  reacted  mainly  to  the  flexural  behavior  of  shell,  are  high  at  the  sections  along 
the  major  and  minor  axes  of  the  elliptical  shell  where  high  bending  moments  created  by  the 
external  pressure  are  experienced.  The  maximum  hoop  stress  Is  87,703  psi  in  compression  for  the 
Inner  fibers  at,the  major  axis  section  and  62,800  psi  in  tension  for  the  outer  fibers.  The  through¬ 
thickness  shear  stress  distribution  is  shown  In  Figure  13.  As  shown,  the  maximum  shear  stress  is 
-6,030  psi,  which  occurs  at  the  neutral  axis  locations  of  the  cross  sections  at  the  transition  area 
of  the  shell  where  D-insert  support  discontinues. 

DISTRIBUTIONS  OF  STRAIN  COMPONENTS 

The  longitudinal  (ej  j),  transverse  (e22).  normal  (€33)  and  through-thickness  shear  (e  13)  strain 
distributions  due  to  the  maximum  pressure  of  459. 1  psi  for  the  upper-bound  property  model  are 
dispatched  in  Figures  14  to  17  respectively.  These  strain  distributions  correspond  to  the  stress 
contour  plots  dispatched  in  Figures  11,  12,  10,  and  13  respectively. 

The  longitudinal  strain  (e  j  j)  distribution  at  a  cross  section  of  the  shell  is  essentially  resulting 
from  the  combined  effects  of  the  flexural  deformation  caused  by  the  pressure  loading  and  the 
stretching  created  by  the  stack  load.  As  shown  In  Figure  14,  the  e*  1  at  the  upper  fibers  of  the 
cross  section  at  the  flat  are  in  the  state  of  high  compression  and  lower  fibers  in  high  tension  due 
to  the  maximum  bending  moment  acting  on  the  section,  while  the  middle-surface  fibers,  not 
affected  by  bending,  arc  in  low  tension  caused  by  the  stack  load.  This  trend  continues  with 
decreasing  flexural  effect,  as  the  bending  moment  decreases  toward  the  transition  point  of  the 
shell.  It  reaches  a  section  in  the  midtransition  zone  where  the  bending  moment  becomes  zero  and 
only  a  small  tensile  strain  caused  by  the  stack  load  remains.  Beyond  this  point,  however,  the 
trend  picks  up  again,  but  with  reverse  effect.  Obviously,  the  bending  moment  has  changed  Us  sign 
at  this  Juncture.  The  reverse  bending  moment  reaches  its  maximum  effect  at  the  end  curve  of  the 


6 


shell  where  the  upper  fibers  of  the  cross  sections  are  in  high  tension  and  innei  libeis  in  hijdi 
compression.  The  effect  of  bending  moment  reversal  can  also  be  observed  by  reviewing  ihc 
defonned  shape  of  the  shell  (Figure  8)--the  flat  of  the  shell  bends  downward  while  the  rum-  bendi- 
outward. 

Since  there  Is  no  direct  loading  effect  in  the  axial  direction  of  the  shell,  the  transverse  strain 
(622)  distribution  shown  in  Figure  15  is  mainly  determined  by  the  Poisson's  effect  of  the 
longitudinal  strains  (e  j  j).  Thus,  the  e22  distribution  is  similar  in  trend,  but  lesser  in  magnitude 
and  opposite  in  direction  to  that  of  j  (Figure  14).  Note  that  a  Poisson's  ratio  of  0.3  was  used  for 
all  three  material  directions  In  the  analysis. 

The  normal  strain  (633)  distribution  is  mainly  determined  by  the  direct  compression  effect 
caused  by  the  pressure  loading,  in  combination  with  the  Poisson's  effect  of  the  longitudinal 
strains.  As  shown  in  Figure  16,  this  results  in  a  normal  strain  distribution  that  Is  predominantly 
compressive  but  retains  some  small  tensile  regions  at  the  upper  fibers  of  the  flat  and  the  inner 
fibers  of  the  curve  where  high  compressive,  longitudinal  strains  caused  by  bending  effects  are  at 
their  maxima.  However,  the  normal  strain  distribution  at  the  inner  surface  region  of  the  curve  is 
greatly  affected  by  the  gapping  created  at  the  shell  and  D-insert  interface  by  the  high  pressure. 

The  gradient  of  the  strain  distribution  at  this  area  is  quite  high,  especially  at  the  area  just  beyond 
the  remaining  D-insert  contact  where  the  maximum  tensile  strain  occurs,  as  shown  in  Figure  16. 

Consistent  with  the  shear  stress  (S 13)  distribution  shown  in  Figure  13.  the  through -thickness 
shear  strain  distribution  (613)  shown  in  Figure  17  indicates  that  the  maximum  shear  strain 
occurs  at  the  neutral  axis  locations  of  the  cross  sections  in  the  transition  area  of  the  shell 

COMPARISON  OF  STRAIN  DATA 

Figure  1  shows  the  locations  and  orientations  of  the  two  types  of  strain  gages  that  were  placed 
on  the  Shell  010  for  strain  measurements  In  the  pressure-plus-stack  load  test.  Gages  of  type  CEA- 
06-062UR-350  are  small  45°  rectangular  single-plane  rosettes  aligned  with  the  circumferential 
and  radial  directions  on  the  cross  sectional  end  of  the  shell  for  longitudinal  (e  j  j),  normal  (633) 
and  through-thickness  shear  (e^  strain  measurements.  Gages  of  type  CEA-13-125WT-350  are 
two-element  90°  "tee"  stacked  rosettes  aligned  with  circumferential  and  axial  directions  on  the 
inside  surface  of  the  shell  for  longitudinal  (ej  j)  and  transverse  {^22)  sbnaln  measurements. 
Eighteen  strain  gages  are  instrumented  at  the  critical  locations  of  the  shell  that  includes  the  flat, 
the  midtransition  and  the  transition  areas,  as  shown  in  Figure  1.  The  calculated  strain  values  at 
the  corresponding  uodeo  of  the  upper  arid  lcwer-’u>.ui  d  property  FEMs  are  compared  with  the 
test  data  reported  by  Maltby  (1991).  As  shown  in  Figures  18  to  35,  the  calculated  strains  for  the 
simulated  pressure  load  Increments  in  load  step  2  are  plotted  against  the  recorded  strain  versus 
stack  stress  curves  for  several  cycles  of  stack  and  pressure  loadings.  Noted  that  the  calculated 
strain  values  for  stack  preloadlng  (load  step  1)  are  omitted  from  the  plots  for  clarity. 

These  figures  show  the  predicted  strain  values  foi  the  upper-bound  property  model  correlated 
better  with  the  test  data  than  the  lower-bound  ones.  Table  3  summarizes  the  quality  of  strain  data 
correlations  between  the  predicted  values  for  the  upper-bound  property  model  and  (he  test 
results. 

In  general,  the  predicted  values  of  strain  components  for  the  upper-bound  property  model 
correlated  very  well  with  the  test  data  at  all  the  locations  considered.  Figures  18  to  35  show  the 
strain  components  calculated  by  FEMs  not  only  have  the  same  signs  as  the  test  data,  but  also 
show  the  similar  trends  with  respect  to  the  increasing  pressure.  However,  the  correlation  for  the 
magnitudes  of  the  interlaminar  shear  strains  (e^)  at  areas  surrounding  the  shell  and  D-insert 
contact  are  mostly  fair  only.  This  discrepancy  may  be  attributed  to  several  reasons  as  follows: 
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1  Uniform  shear  modulus  (O  ^ 3)  was  assumed  throughout  the  entire  thickness  of  the  current 
shell  model.  This  estimated  property  may  be  low  for  the  inner  surface  of  the  composite  shell 
whose  material  properties  may  vary  through  the  thickness 

2.  Figure  17  shows  the  strain  gradient  at  the  inner  shell  surface  area  surrounding  the  I) -insert 
contact  are  steep  due  to  the  strss  concentration  effect  The  rather  coarse  and  uniform  mesh 
of  the  current  FEM  may  not  provide  ABAQUS  with  an  adequate  resolution  to  accurately 
calculate  a  nodal  strain  value,  which  is  averaged  over  the  elements  connected  to  the  node  in 
ABAQUS  analysis. 

3.  Strains  are  calculated  at  a  point  in  the  FEA,  but  they  are  measured  over  a  finite  length 
(0.062  Inch  for  the  rectangular  type  gages)  In  the  test.  An  Inherent  error  between  the  two 
exists. 

Table  3.  Quality  of  Strain  Data  Correlations  Foi  Upper-Bound  Property  Mode! 


Strain 

Node 

Fig. 

ell 

e33 

e13 

e22 

Zone  Gage 

No. 

No. 

(LONG.) 

(NORM.) 

(SHEAR) 

fTRANS.) 

Remarks 

F 


L 

M1IM/M2IM 

1 

18/19 

Fair/ Good 

NA 

NA 

Good 

Biaxial  Gages 

A 

MXAI 

7 

20 

Good 

Very  Good 

Good 

No  shear  strain 

T 

MXAM 

126 

21 

Good 

Good 

Good 

No  shear  strain 

&  others  are  small 

MTXAI 

521 

22 

Good 

Good 

Poor 

MTXAM 

604 

23 

Good 

Very  Good 

Good 

T 

T2XAI/T2XBI 

958 

24/25 

Very  Good 

Very  good 

Good 

R 

T2XAM 

1022 

32 

Good 

Fair 

Fair 

A 

Very 

N 

TD1IM/TD2IM 

963 

26/27 

Very  Good 

NA 

NA 

Good 

Biaxial  Gages 

S 

TDXAI/TDXBI 

969 

28/29 

N A/ Good 

Very  Good 

NA/Falr 

I 

TDXAM/TDXBM 

1026 

30/31 

Very  Good 

Very  Good 

Good/Fair 

T 

I 

O 

TDXAO 

1104 

33 

Fair 

Good 

Good 

T1XAI 

1115 

35 

Good 

Good 

Fair 

N 

T1XAM 

1134 

34 

Very  Good 

Good 

Fair 

Ranking:  Very  Good- -Predicted  values  within  5%  of  Test  Data 
Good--Predicted  values  within  10%  of  Test  Data 
Fair- -Predicted  values  within  25%  of  Test  Data 
Poor— Predicted  values  greaier  than  50%  of  Test  Data 
NA— Test  Data  not  available 

Comments:  (1)  Good  to  very  good  correlations  for  longitudinal  strains  (e  j  j)  at  all  locations, 
except  node  1 104  where  the  correlation  was  only  lair. 

(2)  Good  to  very  good  correlations  for  normal  strains  (033)  at  all  locations, 
except  node  1022  where  the  correlation  was  only  fair. 

(3)  Good  to  fairly  good  correlations  for  shear  strains  (e  j  3)  at  all  locations, 
except  node  521  where  the  correlation  was  poor. 
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CONCLUSIONS  AND  RECOMMENDATIONS 


The  transducer  shell  FEM  constructed  with  ABAQUS'  20-node  solid  elements  using  linear, 
homogeneous,  and  orthot,  ,pic  "smeared"  material  properties  has  predicted  the  global  response 
and  strain  distribution  a*  the  loaded  shell  that  correlate  well  with  the  test  data  obtained  from  the 
stack  load  and  hv  >  jtatlc  pressure  tests.  The  model  should  provide  a  good  base  model  for  the 
sensitivity  stud}  ot  the  shell's  structural  performance  to  thickness  variation  and  to  nonunifomi 
material  property  distribution  in  the  shell  resulting  from  the  variation  in  parameters,  such  as 
moduli  oi  the  constituents,  fiber  volume  fraction  and  ply  arrangement  of  the  laminate.  However, 
the  model  needs  to  be  refined  at  areas  in  contact  with  the  D-insert.  as  well  as  to  simulate 
.•onllnear  material  behaviors  that  were  observed  in  the  pressure  tests.  The  effects  of  variations  in 
longitudinal  and  through-thickness  properties  (E^j  and  G 13)  on  the  flexural  behavior  of  the  shell, 
which  is  vital  to  the  acoustic  performance  of  the  flextensional  transducer,  must  be  investigated 
further. 

As  predicted  by  the  FEMs.  a  gap  formed  at  the  shell  and  D-insert  interface  when  the  shell  was 
pressurized  beyond  400  psi  in  the  tests.  Gapping  has  an  adverse  effect  on  the  acoustic 
performance  of  the  transducer  shell  and  can  reduce  the  fatlque  life  of  the  shell.  It  must  be 
prevented. 

Finally,  a  stress  analysis  cannot  be  concluded  without  mentioning  the  issue  of  failure.  Failure 
criteria  are  more  complex  in  composites,  which  can  delaminate  or  have  either  fiber  or  matrix 
failure.  The  problems  in  defining  failure  modes  for  the  composite  shell  and  the  issue  of  what 
constitutes  a  proper  failure  criterion  in  predicting  reasonable  fatigue  life  for  the  transducer  shell 
must  be  addressed  in  future  work. 
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Figure  1.  Test  configuration  of  Shell  010  showing  strain  gage  locations. 


Uz  =  0  at  Z  =  -2.85-in  Plane 


Figure  2.  Transducer  FEM  with  boundary  and  two-step  loading  conditions 
simulating  stack  load  and  hydrostatic  pressure  tests. 
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Figure  4b.  Second  flexural  mode  of  shell. 


11 


STACK  LOAD  (lb) 


140x10 


SHELL 010 


o  •  - - 


-r—~ — r 


£  legend  ::::: 

^  - TEST  DATA 

±  o  ANALYSS  RESULT  BASED  ON 
"  LOWER  80UN0  MATERIAL  PROPERTIES  — : 

t  A  ANALYSIS  RESULT  BASED  ON 
£  UPPER  BOUND  MATERIAL  PROPERTIES  11" 


MAJOR  AXIS  DISPLACEMENT  (in) 

Figure  5a.  Stack  preload  versus  major  axis  displacement. 
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MINOR  AXIS  DISPLACEMENT  (in) 


Figure  5b  Slack  preload  versus  minor  axis  displacement 


MINOR  AXIS  DISPLACEMENT 
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Figure  9.  Enlarged  view  of  deformed  shell  and  D-insert  interface  under  459-psi  pressure. 
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Figure  12.  Transverse  stress  (S22)  distribution  due  to  459-psi  external  pressure. 
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Figure  13.  Through-thickness  shear  stress  (S,3)  distribution  due  to  459-psi  external  pressure 
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